On quantization of weakly nonlinear lattices. Envelope solitons. 



V. V. Konotoplfl and S. Takenot 
I Departamento de Fisica da Faculdade de Ciencias and Centra de Fisica da Materia Condensada, Universidade de Lisboa, 
Complexo Interdisciplinar, Av. Professor Gama Pinto 2, P-1649-003, Lisbon, Portugal. 
J Faculty of Information Science, Osaka Institute of Technology, 1-79-1 Hirakata, Osaka 573-01, Japan 

(February 1, 2008) 

A way of quantizing weakly nonlinear lattices is proposed. It is based on introducing "pseudo- field" 
operators. In the new formalism quantum envelope solitons together with phonons are regarded as 
elementary quasi-particles making up boson gas. In the classical limit the excitations corresponding 
to frequencies above linear cut-off frequency are reduced to conventional envelope solitons. The 
approach allows one to identify the quantum soliton which is localized in space and understand 
existence of a narrow soliton frequency band. 



PACS number(s): 03.40.Kf 

One of the most important tasks of the modern physics 
of atomic lattices is experimental observation of nonlin- 
ear localized excitations which are beUeved to exist. Al- 
though in specially prepared chains of coupled oscilla- 
tors such excitations have been obtained experimentally 
101 in the microscopic world their observation is hardly 
possible in a direct way. The first report on observa- 
tion of intrinsic localized modes, although in a magnetic 
model, appeared only recently ||^ . This rises a problem of 
construction of statistical mechanics of nonlinear lattices 
which coveres both linear and nonlinear excitations and 
allows one to describe contribution of the localized modes 
to macroscopic (measurable) characteristics of solids. As 
the first step one has to provide quantization of atomic 
chains. 

Anharmonic quantum lattices have been considered in 
literature. Two approaches developed so far one can clas- 
sify following as the number state method ^ and the 
quantum inverse scattering tecnique |4j . Another method 
based on exact numerical diaganolization of a hamilto- 
nian has been developed in |^,p| where breather modes 
were obtained in a lattice with on-site potential |^ and 
in a coupled electron-phonon system |^ , both arc nonin- 
tegrable. The breather modes are identified through the 
spectrum and dynamic correlation functions. 

The mentioned approaches deal with nonlinear quan- 
tum hamiltonians given apriori. An alternative and 
somehow complementary way of introducing quantum 
systems is quantization which starts with a classical 
hamiltonian. As it is evident, the canonical quantization, 
i.e. substitution of c-numbers by their operators satis- 
fying canonical commutation relations, formally can be 
provided in the case of nonlinear lattices, as well. How- 
ever, in a generic situation this way is not tractable nei- 
ther analytically nor numerically. 

At the same time, methods of analytical description 
of nonlinear classical lattices are well elaborated today. 
They are based on introducing different small parame- 
ters and result in various spatially localized excitations. 
Among them we mention envelope solitons (ES) (see e. 



g. ||7|) which appear when the wave amplitude is consid- 
ered as a small parameter and intrinsic localized modes 
or breathers (see e. g. ^). In the case where a small pa- 
rameter can be identified as coupling between neighbor 
sites a quantization can be performed in a way proposed 
in The role of the coupling constant magnitude in 
the context of quantum lattices has been discussed in . 

In the present Letter we concentrate on quantization 
of a nonlinear lattice when the small parameter is the 
nonlinearity, and thus allowing in classical limit existence 
of ES. The respective physical limit can be interpreted 
as opposite to one considered in Also in contrast to 
the previous studies of nonintegrable quantum nonlinear 
lattices, ES appears to be a dynamical object and thus 
allows one to construct a kinetic theory of interacting 
solitons and phonons. 

There are several points to be reflected by the the- 
ory. It has been established that quantum nonlinear lat- 
tices have well pronounced soliton band in the spectrum 
which is extremely narrow for a given nonlinearity 
This means that the respective object could be unam- 
biguously identified either as ES or as intrinsic localized 
mode depending on the value of the frequency. Another 
point is the identification of a quantum ES which must 
go beyond computing the frequency along. Indeed, an 
important difference between ES and a phonon is that 
the soliton is spatially localized while the phonon is lo- 
calized only in the momentum space. Thus we intend to 
obtain a quantum soliton as an object, a quasi-particle, 
localized in space (a possible way of construction of such 
wave packets is described in |^]). This is contrast to the 
quantization used the field theory [l0|] where a kink, be- 
ing a topological object, is considered as a vacuum state. 
Finally, an important property of the classical theory is 
that ES are governed in the leading order by the nonlin- 
ear Schrodinger equation |^ what means that in the lead- 
ing order they do not interact with linear phonons and 
can be regarded as independent quasi-particlcs. Thus in 
the leading order the quantum theory must allow one to 
introduce creation and annihilation operators for quan- 
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turn soliton and neglect its interactions with other quasi- 
particles. In this way the interaction hamiltonian will 
naturally appear as a perturbation and thus the gas of 
interacting quasi-particles - solitons and phonons - can 
be considered within the conventional perturbation tech- 
nique of the field theory. 

It is to be mentioned that the approach developed be- 
low have some points similar to one proposed in Jill for 
quantum optical solitons. In particular, quantum ES in 
PI is described in terms of the field operators (we call 
them pseudo-field operators) which are the Fourier trans- 
form of the conventional creation and annihilation op- 
erators. There are however several essential differences 
between our method and the approach of |pl]] (apart 
from different systems to be quantized). First, we start 
with the original (i.e not approximated) classical dynam- 
ical system. Second, we introduce a soliton as a quasi- 
particle. Third, we consider a highly rarified gas of quasi- 
particles in the quasi-classical limit (while in the 
Hartree approximation for a photon gas was exploited). 
Finally, quantum ES emerges as a matrix element of the 
creation operator applied to the vacuum state JT2| (no- 
tice that in the classical theory envelope solitons are not 
necessarily real). If however one considers an ES as an 
excitation of the displacement field (see w„ below) then 
the quantum ES appears to be an observable. 

Speaking about quantization of ES one can identify 
one essential feature of the quantum model. In the clas- 
sical theory the ES amplitude (and frequency detuning) 
is not a fixed value (see e.g. Q). That is why, the energy 
associated with the classical ES is not a fixed quantity. 
The respective quantum problem has a scale defined by 
the Planck constant: it is natural to expect that the ex- 
citation characterized by the frequency uj will have the 
energy huj. Taking into account that the linear oscillator 
energy is given by moj'^A^ (here to is a mass and A is 
the amplitude) one concludes that the amplitude of the 
quantum ES must be of order of A = ^jTi/iom. On the 
other hand the amplitude of the soliton must be small 
enough. In order to specify this last requirement, let 
us assume that the interaction energy of neighbor sites is 
U{x) and define a characteristic scale of the energy varia- 
tion: L ~ (I j7(jy^^iE^I)^^^^- Then the relation between 
the nonlinear and linear terms in the classical model is 
of order of [A/LY. It is this value that appears to be a 
small parameter in the classical theory, or more precisely 
AjL ^ 1. Then substituting the estimate for the ampli- 
tude obtained above on the basis of quantum approach 
one gets the requirement ^^^2 <C 1 as a condition for the 
validity of the small amplitude expansion. Scaling out all 
the quantities and assuming that the characteristic spa- 
tial size of the excitation is one, L '--^ 1 we conclude that 
the effective small parameter of the problem is ?i, i.e. we 
a re dealing with the quasiclassical limit. 

We consider quantization of a one-dimensional 



monoatomic lattice described by the hamiltonian 
^ = \^^^ \ X! ■^2("-i,n2)a;nia;„2 

(1) 

ni , . . . ,714 

where Xn and p„ — mx„ are displacements of atoms 
from the equilibrium position and their linear momenta, 
K2{n\^ni) and K^^inx^ ■■■,n/^ are linear and nonlinear 
real force constants. Overdot hereafter stands for the 
derivative with respect to time. We concentrate on 
the case of the force coefficients symmetric with re- 
spect to permutation of the indexes, e.g. Kiijii^n^) = 
K2{n2,ni). Also it is assumed that J2ni ^2("'i,«2) — 0. 

First of all we recall some facts of the theory of linear 
lattices. Consider vibrations of a one-dimensional (pos- 
sibly disordered) lattice described by the Hamiltonian 

n " n,l 

where to„ is the mass of the nth atom and K{n, I) = 
K(l, n) are real force constants. It is convenient to intro- 
duce new dependent variables ii„ = ^/m^Xn and force 
constants J{l,n) — K{1, n) / ^to;to„. The respective 
equation of motion reads 

Un+^J{n,l)ui = Q (3) 

It will be assumed that the lattice consists of A/", 
A/" >• 1, atoms and is subject to the cyclic boundary 
conditions. Then equation (H) can be associated with 
the linear spectral problem [^q{n) — ipq^n + Af)] 

Y,J{nMS)=^l%{n), (4) 
I 

where real eigenvalues are squared eigenfrequencies 
and q denotes eigenmodes, eigenfunctions ipq(ri) and 
i/jg (n) correspond to the same eigenfrequency (an overbar 
stands for complex conjugation). Introducing the matrix 
J with the element J{n, I) placed at nth row and Ith col- 
umn, a ket-vector \q} = col{.. .,ipq{n — l),'ipq{n),ipq{n + 
1), ...) and a bra-vector {q\ — (..., ipq^n—l), 'ipq(n),ipq{n + 
1),...) we rewrite the eigenvalue problem (^ in the 
form 3\q} = ujqlq}- We assume that ujq — lu-q, i. 
e. i^q{n) = ip-q{n). ipq(n) constitute an orthonor- 
mal complete set {gig'} = (") ~ ^9-9'' ^^'^ 

Tliqi'q{''^)'^q{^) ~ ^n,i (here 5i^n is the Kronecker delta). 

Quantization of lattice (|^) can be reduced, first, to 
the problem of the diagonalization of the matrix J 
and, second, to the quantization of independent lin- 
ear oscillators. In the normalized variables the ki- 
netic energy is written in the form {7r|I|7r} where 
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\n} =col(..., 7r„_i, 7r„, 7r„+i, ...), 7r„ = is the mo- 
mentum conjugate to u„, and I is the unit ma- 
trix. Thus any matrix diagonahzing the potential 
energy will preserve the diagonal form of the ki- 
netic energy. To provide diagonalization of the po- 
tential energy we construct a A/" x TV matrix 'if = 
(..., \qj-i}, \qj}, \qj+i}, ...) [here numbering of the discrete 
eigenvalues is introduced]. As is clear J'S' = '3/0^, where 
17 =diag{...,LUq^_-^,ujg.,ujq.^-^, ...), i.e. the matrix * diag- 
onalizes J: 'd'^^J* = fi^. If one represents \u} = 
and {u\ = where \u} =col(..., m„, Un+i, ■■■) 

and \X} =co\{..., Xq._-^, Xq. , Xq.^-^, ...) is a column ma- 
trix of the normal coordinates with Xq = X^q, the po- 
tential energy will take the form {X|r2^|X}. Finally 



H^l{p\p} + l{x\n'\x}., 



(5) 



where \P} ~co\{..., Pq.-^, Pq. , Pq.^-^, ...) is related to the 
representation |7r} = 'S'l^'}- 

In order to quantize the linear disordered lattice one 
introduces phonon annihilation and creation operators 



n - r 'if 



\/2hujq 

i 



Pn 



(6) 



where Xq and Pq are operators of the c-numbers Xq and 
Pq satisfying the canonical commutation relations. Then 
the Hamiltonian is reduced to the conventional form Hi = 
hojq {ci'qdq + i) . Opcrators of displacements and of 
linear momenta read 



[i)q{n)aq + ^:q{n)a\] 




-[i>q{n)aq - i/'g("-)a, 



(7) 



Let us introduce Schrodinger pseudo-field operators ac- 
cording to the relations 



(8) 



They describe creation and annihilation of a quasi- 
particle at the site n. We emphasize the main distin- 
guishing feature of these operators compared with the 
canonical field operators: '4^q{n) is an eigenfunction of 
the classical problem rather than a wave-function. 
It is not difficult to verify that 



n n 

and there exist commutation relations 



(9) 



The operator of number of quasi-particles is given by 

^ = E«9^? = E^"V'n (10) 
q n 

Finally, following the conventional procedure we intro- 
duce Heisenberg pseudo-field operator 



i , ^ Ht\ - Ht 

w„(i) = exp I I— I ipn exp I 



It solves the equation 



& 



(11) 



(12) 



Let us now consider a nonlinear lattice with classical 
Hamiltonian (^) after the renormalization J2irii,n2) = 
K2{ni, 712 )/m and J4(rii, n^) = K4{ni, n/^jvr? . To 
obtain explicit form of the quantum hamiltonian we as- 
sume that the definition holds in the nonlinear case. 
We substitute operators of the displacement and linear 
momenta into the Hamiltonian, and express the re- 
sult through the pseudo-field operators, rearranging the 
last ones in the normal order and dropping the constant 
which corresponds to the energy of the lattice vacuum. 
Then Tti — Hq + Hint ■ The operator H^ has the form 

-ffo = E '5'ni,«2*ni*n2 = E ^^^^ 



with the kernel Sni,n2 — ^'^g''/'g("'i)'0g("2)- Now 
ipq^n) is an eigenfunction and ujq is a frequency of the 
"nonlinear" eigenvalue problem. The " nonlinearity con- 
tribution" emerges from the ordering procedure. More 
precisely, the eigenvalue problem (^) is now considered 
with the kernel 



J{n, I) = J2(n, I) + Jd{n, I) 



(14) 



where J2(n, I) is the force constant of the underline linear 
lattice and the "deformation" Jd{n, I) is given by 



Jd{ni,n2) 



3h 



E E '^4("-l,«2,^l,^2)'0g(^l)V'9('2) 



9 hj2 



The interaction potential Hint has the form 

n 

+ X&t a, x| _1_ \|t \T,t 1 it fl;t »T;t 1 

' *n2 I *Tii *n2 ^"4 ' ^ rii * Ti2 ^ 113 ^ n4 J 



with the kernel 



Sn-ii 



Ji{li, I2, h, h) 

16 V^«lWg,W,3tJq, 



^qi{h)-il^q2{h) X 



i'qi {h)i^qi {k)'ipqi ("-l)'092 (»^2)'093 ("-3)'(/'g4 (^^4) 
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In the small amplitude limit which corresponds to 
small J4, Hint can be considered as a perturbation. That 
is why we introduce pseudo-field operators in the inter- 
action representation, '^n\t), by formula ( |ll]) where Ha 
is used instead of H. Then one can ver ify that ¥n\t) 
solves (|l|) with J{n,l) given by (|lj). 

In order to show how ES comes out from the above 
approach we notice that the ES solution corresponds to 
the situation when only one quasi-particle is excited in 
the lattice. In what follows quasi-particles corresponding 
to the frequencies inside the spectrum band of the lin- 
ear lattice, Lo e [0,wo) where ujq is the cut-off frequency 
of the underline linear lattice, are called phonons while 
excitations with lo > ujo will be referred to as solitons. 
Thenl' = \E'p,, + *s- 

Let us introduce the notation jn^j^ , for the wave 
function of a state when nq^ quasi-particles with the wave 
number qj are excited. Then the "one soliton" state, 
when there exists only one eigenvalue ujs bigger than luq, 
is |0, 0; 1) and ipsin) = (0, 0; 0|^|0, 0; 1) [we use 
the notation ipsin) = '0q=Wa, i^s = ^q=TT/a]- In order to 



obtain 1}}^ we rewrite Eq. in the form 

iolil)s{n) = ^ iljs{m)J2{n, m) + — X 

J4(n, mi,TO2,m3)V's(TOi)?/is(m2)?/'s("l3) (15) 

mi ,7712 ,m3 

As far as the small amplitude limit is under consideration 
the frequency LUg is close to the cut-off frequency of the 
underline linear lattice, loq, and one can introduce a small 
parameter e (e ^ 1) through the relation cj^ — cjq = 
e^u>Q where = 2{uJs/ll!q — 1) is the frequency detuning 
towards the forbidden zone. 

We look for the solution of ( [l5| ) in the form 



ipsin) 



3n 



(16) 



where (fn is the eigenvalue of the problem for the perfect 
linear lattice, i.e. J{n,l) = J2{n,l) in (^. Then expand- 
ing A{e{n — I)) (here I depends on the number of nearest 
neighbors interacting with given atom) in the Taylor se- 
ries with respect to e/, introducing x = tan (treated as a 
continuum spatial variable) and taking into account that 
Ipsin — I) — Ipsin + I) one arrives at 



J4(ni, ...,ni)LPmVn2Vn3VnAA\'^ A = 



which is obtained for monochromatic solution of the ID 
nonlinear monoatomic lattice within the framework of 
the conventional multiscale analysis 0. 



The obtained results corraborate with the known ones. 
In particular, it has been shown that for the given non- 
linearity the frequency of the quantum soliton is fixed by 
the normalization condition [ pl| . This condition defines 
e introduced in (|l^) and must be computed for each lat- 
tice separately. It is however clear that in any case the 
result gives e ~ ^/h which means the consistency of the 
expansion. On the other hand this result does not con- 
tradict to the existence of a narrow soliton band which 
has been explained in |^,^,|ll| with help of different argu- 
ments. Indeed, quantum ES is characterized not only by 
the "carrier wave" frequency but also by the frequency 



smearing which is of order of (?i/m) 



. Our re- 



sults corraborate also with recent findings in the classi- 
cal statistical mechanics of nonlinear systems. Namely 
it is well known that boson display a tendency of cre- 
ating clusters (see e.g. Q). This means that creation 
of a phonon and a quantum ES sufhciently close to each 
other will end up in clustering of the two quasi-particles. 
This can be viewed as absorption of the phonon by the 
soliton. This phenomenon has been observed recently in 
numerical experiments . Finally, in the leading order 
with respect to the nonlinearity solitons and phonons are 
noninteracting objects. Interaction is taken into account 
in higher orders of the expansion. 
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